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ABSTRACT 

We generally analyze prompt high-energy emission above a few hundreds of GeV due to synchrotron 
self-Compton scattering in internal shocks. However, such photons cannot be detected because they 
may collide with cosmic infrared background photons, leading to electron/positron pair production. 
Inverse-Compton scattering of the resulting electron/positron pairs off cosmic microwave background 
photons will produce delayed Me V- GeV emission, which may be much stronger than a typical high- 
energy afterglow in the external shock model. We expand on the Cheng & Cheng model by deriving the 
emission spectrum and duration in the standard fireball shock model. A typical duration of the emission 
is ~ 10 3 seconds, and the time-integrated scattered photon spectrum is ^~(p+ 6 )/ 4 , where p is the index 
of the electron energy distribution behind internal shocks. This is slightly harder than the synchrotron 
photon spectrum, zy _ (?>+ 2 )/ 2 . The lower energy property of the scattered photon spectrum is dependent 
on the spectral energy distribution of the cosmic infrared background radiation. Therefore, future 
observations on such delayed Me V- GeV emission and the higher-energy spectral cutoff by the Gamma- 
Ray Large Area Space Telescope (GLAST) would provide a probe of the cosmic infrared background 
radiation. 

Subject headings: diffuse radiation — gamma rays: bursts — gamma rays: theory 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are the brightest electro- 
magnetic phenomena in the universe. Their radiation fea- 
tures have been well understood even though their origin 
has been unknown. In the standard fireball shock model 
(Piran 1999; van Paradijs, Kouveliotou & Wijers 2000; 
Meszaros 2002), GRBs are explained to be due to the dis- 
sipation of kinetic energy of an expanding fireball with 
an average Lorentz factor more than 10 2 in the internal 
shocks produced by collisions between different shells in 
the fireball. Afterglows are also considered to be due to 
the dissipation of kinetic energy in the external shocks 
generated by collision of the fireball with its surrounding 
medium. Such a model has been supported by a variety of 
observations on GRBs and afterglows (e.g., Wijers, Rees & 
Meszaros 1997; Waxman 1997; Vietri 1997). Theoretically, 
it has also predicted high-energy emission from GRBs and 
afterglows (for a review see Meszaros 2002). The observed 
spectra of GRBs sometimes extend to GeV-TeV photon 
energies. For example, the EGRET experiment detected 
an 18 GeV photon from GRB 940217 (Hurley et al. 1994), 
and the Milagrito air shower experiment tentatively de- 
tected the TeV emission from GRB 970417a (Atkins et al. 
2000). These GRBs might have originated at a redshift of 
z < 2 (Salamon & Stecker 1998; Totani 2000). 

The purpose of this paper is to present a new prediction 
of delayed MeV-GeV emission of GRBs in the cosmic back- 
ground radiation fields. Such emission arises from inverse- 
Compton scattering of electron/positron pairs off cosmic 
microwave background (CMB) photons. In §2 we generally 
discuss prompt high-energy emission above a few hundreds 
of GeV due to synchrotron self-Compton scattering in in- 
ternal shocks based on the standard fireball model. How- 



ever, such high-energy photons cannot be detected because 
they may be absorbed in interacting with the cosmic in- 
frared background radiation (Stecker, De Jager & Salamon 
1992; Madau & Phinney 1996; MacMinn & Primack 1996; 
Malkan & Stecker 1998; Salamon & Stecker 1998). In §3 
we show that inverse- Compton scattering of the resulting 
electron/positron pairs off CMB photons will produce de- 
layed MeV-GeV emission by deriving the time-integrated 
spectrum of such emission and by estimating the emission 
duration. In §4 we summarize our findings and discuss 
their implications. Cheng & Cheng (1996) have discussed 
similar delayed MeV-GeV emission, and here we further 
calculate the emission spectrum and estimate the emission 
duration in the standard fireball shock model. 



2. HIGH-ENERGY RADIATION FROM INTERNAL SHOCKS 

Let us consider a compact source that produces a wind 
(i.e., a variable fireball) with an average luminosity of L 
and with an average mass loss rate of M = L/rjc 2 . Ini- 
tially, the bulk Lorentz factor of the wind, T, increases 
linearly with radius, until most of the wind energy is con- 
verted to kinetic energy and T saturates at T ~ r\. Af- 
ter this time, fluctuations of T due to variabilities in L 
and/or M on time scale At could lead to internal shocks 
in the wind at radius R ~ T 2 cAt. We also assume that 
the internal shock Lorentz factors Ti ~ a few in the wind 
rest frame, so that the internal shocks could reconvert a 
significant fraction of the kinetic energy to internal en- 
ergy, which is then emitted as 7-rays by synchrotron and 
inverse-Compton (IC) radiation of the electrons acceler- 
ated by the shocks. 

The comoving electron number density of the unshocked 
wind matter at radius R is n' e — L/(4irr 6 At 2 m p c 5 ), 
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where m p is the proton mass. The internal energy den- 
sity of the shocked wind matter is thus approximated 
by e' w 4r,(r, - l)n' e m p c 2 (Blandford & McKee 1976). 
The electron energy distribution behind the shocks is usu- 
ally a power law: due/d'ye oc 7~ p for "f e > j m . As- 
suming that e e and cb are constant fractions of the in- 
ternal energy density going into the electrons and the 
magnetic field respectively, we obtain the minimum elec- 
tron Lorentz factor j m — \{p — 2)/{p — l)](m p /m e )e e Ti 
(where m e is the electron mass) and the magnetic field 
strength of the shocked wind matter B = (S^ese') 1 / 2 = 
6.2 x 10 3 e 1 J 2 _ 2 (T l /2) 1 / 2 r^L 1 5 / 2 2 AtZ 1 2 G, where e Bj _ 2 = 
es/10- 2 , r 500 = T/500, L 52 = L/10 52 crgs-\ and At_ 2 = 
Ai/10~ 2 s. Thus, we estimate the characteristic frequency 
of synchrotron radiation from the shock-accelerated elec- 
trons, 

= 3.0xl0 19 [(p-2)/(p-l)] 2 

2 1/2 
X e e,0.5 £ B,-2^ 



(Tj2f/^ lL l 5 / 2 2 AtZ 1 2 li z , (1) 



where e ei0 .5 = e e /0.5. According to Sari, Piran & Narayan 
(1998), the cooling Lorentz factor j c <C j m for typical pa- 
rameters of the internal shock model, implying that the 
shock-accelerated electrons are in the fast cooling regime. 
Therefore, the synchrotron spectral luminosity is 
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At 



Lv m {vlv m ) 1/2 if v < is„ 
L Vm {v/v m )-P/ 2 \iv>v. n 



(2) 



where h is Planck's constant and N v is the synchrotron 
photon number per unit frequency in At. 

Although GRBs are optically thin to electron scattering, 
some synchrotron photons will Compton scatter on the 
shock-accelerated electrons, producing an additional IC 
component at higher energies (Papathanassiou & Meszaros 
1996; Pilla & Loeb 1998; Panaitescu & Meszaros 2000). 
The ratio of the IC to synchrotron radiation luminosity 
can be estimated by 
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where we have assumed that e e 3> £_b and the radiating 
electrons are in the fast cooling regime (Panaitescu & Ku- 
mar 2000; Sari & Esin 2001; Zhang & Meszaros 2001). 
Since the characteristic IC frequency 
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the ratio of the IC spectral luminosity at to syn- 
chrotron spectral luminosity at v m is 



L}? 



7, 
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which is derived from equation (3), and the ratio of the IC 
to synchrotron photon number per unit frequency in At is 



further given by 

_ _ 1/ 



• ^ 2 Y, for v>v, 
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The intrinsic cutoff energy in the IC component is de- 
fined by £™ t = mm(£^ N ,£^), where is the Klcin- 
Nishina limit and £ n is the energy at which a photon 
may be attenuated due to pair production through inter- 
actions with softer photons (also radiated from the inter- 
nal shocks) whose frequency is equal to or larger than 
u an = (Tm e c 2 ) 2 /(h£^). We first estimate £^ N . There 
is the maximum Lorentz factor 7m in the electron energy 
distribution behind the shocks due to the fact that the 
shock-acceleration time cannot exceed the cooling time. 
This gives jm — {3e/[orB(l + V)]} 1 / 2 where e is the elec- 
tron charge and <jt is the Thomson cross section. Thus, 
we obtain the Klcin-Nishina limit £j^ = ^M^nieC 2 ~ 
56e e :o / 5 4 (r 4 /2)- 1 /4 r 5/2 L -i/4 At V2 TeV ^ Tq computc gin 

we adopt an analytical approach proposed by Lithwick & 
Sari (2001). Integrating the synchrotron spectral luminos- 
ity over frequency should lead to the total synchrotron 
luminosity, which further gives the synchrotron photon 
number per unit frequency at v m : N Vm = {(p — 2)/[2(p — 

1) ]}e e (l + Y)- 1 (hv 2 m )- 1 LAt. The total synchrotron pho- 
ton number at frequencies larger than v can be estimated 
by N >v ~ C N » dv = C N Vm { V / Vm )-^+ 2 )l 2 dv = {(p - 

2) /[p(p-l)]}e e (l + Y)- 1 (hv m )- 1 LAt(v/v m )-vl 2 . The at- 
tenuation optical depth is 
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The requirement that t™ = 1 leads to the upper energy 
limit 



= 160e; i2 r 3)/p e B 1/2 2 (r t /2)- 5 ^/^ 



-B,-2\ 
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4(p+l)/p r -(p+2)/(2p) 
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At_ 2 GeV, 
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where the numerical value is quoted for p = 2.2. It is easy 
to see that £g^ ^> £™ for typical parameters and thus the 
intrinsic cutoff energy is in fact given by equation (8). 

3. PAIR PRODUCTION AND RERADIATION IN COSMIC 
BACKGROUND PHOTON FIELDS 

High-energy 7-rays emitted from the internal shocks 
may not only be intrinsically attenuated due to pair pro- 
duction through interactions with softer photons from the 
same radiation regions but may also be absorbed in the ex- 
ternal background radiation fields when these 7-rays travel 
towards the observer. In the latter case, the observed 
cutoff energy, £°^ t , is determined by the redshift (z) and 
comoving energy, £j(z), of the background photons that 
dominate the pair production optical depth (r™) because 
the cross section peaks when £°^ t £ 7 (z) ~ (m e c 2 ) 2 /(l + z). 
It is general that at low redshifts the spectral energy dis- 
tribution of the starlight background radiation peaks at 
infrared wavelengths and thus £°^ t is in the TeV energy 
range (Salamon & Stecker 1998). Detailed predictions of 
£°^t differ from model to model, depending on the theo- 
retical treatment adopted for the stellar emissivity. For 
example, the observed cutoff energy when = 1 for 
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z = 1 is 50 GeV £ £° b t £ 80 GeV for different models used 
in Salamon & Stecker (1998). 

The pair production optical depth, r°*, strongly de- 
pends on the 7-ray energy (£° b ). Salamon & Stecker 
(1998) numerically calculated as a function of the pho- 
ton energy for several fixed redshifts by considering the 
stellar emissivity with and without metallicity correction, 
shown in their Figures 6 and 7, respectively. It can be 
seen from their Figure 6 that if z = 1, then ~ 1 for 
E° h = 50 GeV, but ~ 10 for £° b = 300 GeV. In this 
model of Salamon & Stecker (1998), therefore, a photon 
with energy of £° h ;> 300 GeV must have a pair production 
optical depth » 1. This implies that such high-energy 
7-rays may be locally attenuated once they are radiated 
from internal shocks. For simplicity, we will neglect any 
rcdshift correction for low-rcdshift sources discussed in this 
paper. 

The resulting electron/positron pairs have Lorentz fac- 
tors of 7e EE £° b /(2m e c 2 ) > 7min ee £^/{2m e c 2 ) = 
3 x 10 5 £q X 3oo ! where £§ x = £ § x 300 x 300 GeV is the minimum 
energy of photons that are locally attenuated in the ex- 
ternal infrared background radiation field. The pairs will 
Compton scatter the CMB photons. As a result, the initial 
energy of a microwave photon, Iivq, is boosted by IC scat- 
tering up to an average value ~ 7e^o > 57(£q* 300 ) 2 MeV, 
where hvo — 2.7kT is the mean energy of the CMB pho- 
tons with T ~ 2.73 K and k is the Boltzmann constant. 
The IC lifetime (in the local rest frame) of an electron 
with Lorentz factor of 7e reads 



T(7e) 



3m P c 



47 e (7TUcmb 



< 2.4 x 10 14 (f ™30o) _1 s, (9) 



where w cm b = aT 4 is the CMB energy density and a is 
the radiation constant. It is not difficult to find that the 
typical length in which most of the electron energy is lost, 
~ cr( 7e ), is much less than the distance from the source to 
the observer, implying that energy loss of the electron due 
to IC scattering is also local. Beyond this length, energy 
loss of the electron and its emission become insignificant 
in the absence of any acceleration. 

The typical duration of the scattered photons after 
the GRB trigger in the observer's frame is estimated by 
r ob = max(-r° b , r^), where r° b is the observed IC cooling 
lifetime, 



r ob ~^)< L3x 10 3 (£ - 300 



(10) 



^ob 



^pair/(27gc) is the angular timescale (Pi- 

~ 5.8 x 



and 

ran 1999). Here i? pair = (0.26cr T n m ) 



10 24 (niR,/l cm -3 ) -1 cm is the typical pair-production ra- 
dius (where Uir ~ 1 cm -3 is the cosmic infrared photon 
number density, Protheroe & Stanev 1993). Thus, we have 

r 2 ob < 1.0 x 10 3 (n IR /l cm" 3 )- 1 ^)" 2 s. (11) 

Therefore, most of the scattered photons will reach the 
observer in ~ 10 3 seconds. This time is much longer than 
a typical GRB duration. 

We next derive the scattered photon spectrum. From 
equations (2), (3) and (6), we first write the spectrum of 



-(p+2)/2 



the electrons and positrons as follows 

dN * 9 P-2 V AT ( hv m V /2 . 

if 7min < 7 e < 7max, (12) 

where 7ma x = £o n /(2T7i e c ). The scattered photon number 
per unit frequency is given by (Blumcnthal & Gould 1970, 
hereafter BG70) 



dN~, 



dt'dv 



T( 7e )ded7e 



(13) 



where dN let€ / dt'dv, expressed by equation (2.42) of BG70 
in the Thomson limit, is the spectrum of photons scattered 
by an electron with Lorentz factor of 7e from a segment of 
the CMB photon gas of differential number density n(e)de, 
and t 1 is the time measured in the local rest frame. Please 
note that the scattered photon number per unit frequency 
given by equation (13) is time- integrated by multiplying 
r( 7e ) in equation (2.61) of BG70. After integrating over -j e 
and e, we obtain the ratio of the scattered to synchrotron 
photon number per unit frequency 



Nl 



r7; 



p-2y 



9(m e c 2 )-(P- 2 >/ 2 

2 (p+8)/2 7r 2( ?k .)3 ar 4 

x F( g )(fcT)(P+ 14 )/ 4 (/u/)(P" 2 )/ 4 



2 -p/2y 



F(q) 



7 i(fcT)(H 



"1 (p-2)/4 



(m e c 2 ) 2 
= (p + 4)/2, 



(14) 



F(3) 

where F(q) is a function of q 

F <" " (7^T^ r[(, ' +5,/2KI(5+5)/21 (15) 

with the Gamma function and Riemann zeta function. 
Therefore, the scattered photon spectrum becomes 



N sc K zy -( P +6)/4_ 

The lower energy limit of the scattered photons is 

hv x ~ 7 ^ in (2.7fcT) = 57(f cx 300 ) 2 MeV, 
while the upper energy limit is 

^2^7max(2.7fcT) 

c 16eJ.f - 3)/P ^:_ 2 (r i /2)- 5 ^- 2 )^ 

"(P+2)/p , 



(16) 



(17) 



r 8(p+l)/pr- 

1 500 -^52 



"At* 2 MeV, 



(18) 

which is derived from equation (8). 

We now consider a simple example. The typical pa- 
rameters of the internal shock model are taken: e e> o,5 = 
cb,-2 = Ui = At_ 2 = r</2 = 1, r 500 = 1.4 and p = 2.2. 
Furthermore, £q* 30 o = 1 is assumed. Therefore, we obtain 
hv\ ~ 57 MeV and hv 2 ~ 0.8 GeV. From equation (14), we 
further find iV^ c /N v increases from ~ 5.0 to ~ 5.8 when the 
observed energy of the scattered photons increases from 
hv\ to hv 2 . This leads to a delayed MeV-GeV emission 
component. It is interesting to note that the total energy 
release of such emission is not only comparable to that of 
a typical GRB, but also much larger than the total MeV- 
GeV energy release in an early afterglow as calculated by 
Dcrmer, Chiang & Mitman (2000) and Zhang & Meszaros 
(2001). 
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4. DISCUSSION AND CONCLUSIONS 

We have discussed prompt high-energy emission from 
the internal shocks, which are produced by collisions be- 
tween shells with different Lorentz factors in an ultrarel- 
ativistic wind ejected from the central engine of GRBs. 
Such emission may arise from synchrotron self-Compton 
scattering in the internal shocks. We derived the intrinsic 
cutoff energy due to pair production through interactions 
of high-energy photons with softer photons from the same 
radiation regions. Our intrinsic cutoff energy is consistent 
with the observed spectra of low-redshift GRBs such as 
GRB 940217 and GRB 970417a. However, another cut- 
off energy appears in the 7-ray spectra. This cutoff results 
from interactions of high-energy photons with external cos- 
mic infrared background photons. The electron/positron 
pairs produced during such interactions own a significant 
fraction of the explosion energy and they will Compton 
scatter the CMB photons. We expanded on the Cheng 
& Cheng (1996) model by deriving the emission spectrum 
and duration in the standard fireball shock model. The re- 
sulting emission will be able to reach the observer in ~ 10 3 
seconds. The time-integrated scattered photon spectrum 
is A^ c oc j/~(p+ 6 )/ 4 . This is slightly harder than the inter- 
nal shock emission spectrum, N v oc ^~(p+ 2 )/ 2 . For typical 
parameters of the internal shock model, the observed en- 
ergies of the scattered photons are in the MeV-GeV range, 
and the ratio of the scattered to synchrotron photon num- 
ber per unit frequency, N^ c /N Vl is typically a few. 

Of course, whether or not there is such delayed emission 
depends on the parameters of the internal shock model. 
Some parameters (e.g., L, e e and cb) adopted in this paper 
are consistent with detailed fits to the multi-wavelength 
data of a few afterglows (Wijers & Galama 1999; Freedman 
& Waxman 2001; Panaitescu & Kumar 2002). Two other 
parameters (e.g., p and At) are consistent with the BATSE 
observations (Preece et al. 2000). The bulk Lorentz fac- 
tor of the wind, T, is a crucial parameter because 5™ is 
strongly dependent of T (see equation [8] ) . The observed 
TeV emission from GRB 970417a requires that T be equal 
to or larger than 700, provided that the other parameters 
are given in the simple example o/§5. Furthermore, some 
arguments on the optical flash of GRB 990123 in the re- 
verse shock model (Wang, Dai & Lu 2000; Soderberg & 
Ramirez-Ruiz 2002) show that V may be of the order of 



10 3 (but also see Sari & Piran [1999], who suggested that 
P in this burst could be about 200). If T is close to or 
larger than 600, we can indeed see a delayed MeV-GeV 
emission spectrum. If T < 600, however, the intrinsic cut- 
off energy could be less than the external pair-production 
minimum energy at which a photon is locally attenuated 
through interaction with the infrared background radia- 
tion, so that there could not be the delayed MeV-GeV 
emission discussed here. 

It is seen from equations (6) and (14) that the Comp- 
ton parameter Y determines the TeV emission from syn- 
chrotron self-Compton scattering in internal shocks, and 
thus the spectral photon number of the delayed MeV- 
GeV emission. In addition, additional TeV photons 
produced by other mechanisms, e.g., photo-pion and 
inelastic proton-neutron collisions (Waxman & Bahcall 
1997; Bottcher & Dermer 1998; Derishev, Kocharovsky 
& Kocharovsky 1999; Bahcall & Meszaros 2000; Dai & Lu 
2001), and synchrotron self-Compton scattering both in 
external reverse shocks (Wang, Dai & Lu 2001a, 2001b) 
and in early-time forward shocks (Dermer, Chiang & 
Mitman 2000; Zhang & Meszaros 2001), may be ab- 
sorbed in the cosmic infrared photon field, leading to elec- 
tron/positron pairs. The inverse Compton scattering of 
such pairs off CMB photons may have a non-negligible 
contribution to the delayed MeV-GeV emission. 

The MeV-GeV emission studied in this paper will be 
detectable by next-generation 7-ray satellites such as the 
Gamma-Ray Large Area Space Telescope (GLAST). In the 
Swift-GLAST era, a plenty of GRBs localized by the Swift 
satellite are expected to be detected by GLAST. Such ob- 
servations, particularly on the delayed MeV-GeV emission 
and higher-energy spectral cutoff, may provide a probe of 
the cosmic infrared background radiation. This in turn 
may help to constrain some of the most fundamental un- 
certainties in physical models of the star formation. 
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